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1.  Introduction 


1.A  Background: 

Each  year  about  1 86,000  women  contract  breast  cancer,  of  whom 
46,000  will  die  of  it.  Breast  cancer  is  the  leading  cause  of  cancer 
deaths  amongst  American  women  aged  35  to  50  and  also  has  the 
highest  cancer  mortality  rate  for  younger  women.  Mammography 
is  most  widely  used  modality  for  early  detection  of  breast  cancer. 
The  detection  of  tiny  carcinomas  and  low  contrast  tissue  is  essen¬ 
tial  for  mammography  to  reduce  the  mortality  rate1, 2.  It  has  been 
shown  that  about  20%  of  detectable  cancers  are  overlooked  or  mis¬ 
diagnosed  on  first  detection.  We  obviously  need  a  major  improve¬ 
ment  in  the  quality  of  mammograghic  imaging.  X-Ray  capillary 
optics  has  this  ability. 


Figure  1  Complete  optic  con¬ 
structed  by  configuring  fibers 
through  metal  grids. 


1.B  Benefits  of  digital  mammography 


Conventional  film  screen  mammography  suffers  from  limited  dynamic 
range  and  film  granularity,  which  can  reduce  the  sensitivity  to  detect  micro¬ 
calcifications.  The  primary  theoretical  limitations  of  mammography  are  the 
system  resolution,  which  determines  the  minimum  size  of  the  detectable 
malignancy,  and  the  need  to  expose  the  patient  to  ionizing  radiation.  Using 
innovative  new  technology  --  polycapillary  optics  —  to  improve  system 
resolution  and  reduce  required  dose  will  increase  the  effectiveness  of  this 
proven  screening  modality,  with  a  direct  and  immediate  impact  on  mortality. 
By  matching  polycapillary  optics  to  digital  detectors,  we  can  develop  direct 
digital  mammographic  systems  which  can  avail  themselves  of  the  advan¬ 
tages  of  digital  processing,  including  improved  image  contrast  and  resolu¬ 
tion  at  reduced  radiation  dose.  In  practice,  mammographic  imaging  is  often 
limited  by  quality  assurance  issues,  which  can  also  be  favorably  addressed 
by  digital  processing.  Digital  detection  can  provide  high  dynamic  range 
which,  in  addition  to  improving  contrast,  greatly  increases  the  tolerance  of  the  final  image  to  under  or 
over  exposure.  Digital  images  can  be  enhanced  and  are  amenable  to  computer-aided  diagnosis.  Finally, 
digital  images  can  also  be  transported  quickly  for  skilled  consultations. 


Figure  2  Cross  section 
of  polycapillary  fiber. 
0. 5  mm  in  diameter. 
SEM  photo. 


1. C  What  capillary  x-ray  optics  can  do 

Capillary  optics,  consisting  of  arrays  of  hollow  glass  tubes,  is  a  relatively  new  technology  for  controlling 
X-ray  beams.  X-rays  incident  on  the  interior  of  the  glass  tubes  at  small  angles  are  guided  down  the  tubes 
by  total  external  reflection.  Arrays  of  curved/tapered  capillaries  can  be  used  to  focus,  collimate  and  filter 
X-ray  radiation3,4,5,6’7.  Such  arrays  are  manufactured  by  stringing  hollow  glass  polycapillaries  through 
metal  grids  (as  in  Figure  1)  or  as  a  monolithic  optic. 

The  use  of  Kumakhov  capillary  optics  in  place  of  conventional  scatter  reduction  grids  in  a  mam¬ 
mographic  system  has  significant  potential  to  provide  improved  resolution,  increased  contrast  enhance- 
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ment  and  reduced  dose  in  mammographic  imaging.  The  optics  can  also  be  used  to  mate  the  radiographic 
image  with  a  digital  detector  by  appropriate  choice  of  demagnification  and  separation  to  discrete  chips. 
In  addition,  a  pre-patient  optic  could  be  employed  to  increase  the  available  intensity  in  a  fan  beam  rela¬ 
tive  to  simple  slot  collimating. 


The  fibers  in  this  research  varied  from  0.3  mm  to  4  mm  in  outer  diameter,  with  channel  sizes  (the 
diameters  of  the  holes  within  the  fibers)  ranging  from  4  to  22  pm,  as  shown  in  Table  1  in  section  2.E. 
Open  area  in  the  table  refers  to  the  fraction  of  the  front  face  of  the  fiber  which  is  open  channels,  as 
opposed  to  glass  walls.  No  fiber  can  have  a  transmission  that  is  greater  than  the  open  area. 


X-ray  beams  can  be  bent,  focused  or  collimated  by  a  carefully  curved  polycapillary  array.  X-ray  photon 
energies  are  much  larger  than  the  plasma  frequencies  of  glasses,  which  are  tens  of  electron  volts.  In  this 
regime,  the  real  part  of  the  index  of  refraction  of  glass  can  be  simply  approximated  by 


where  n  is  the  index  of  refraction,  e  is  the  dielectric  constant  of  the  glass,  co  is  the  photon  frequency,  cop 
is  the  plasma  frequency  of  the  material,  and  so  is  the  dielectric  constant  for  vacuum. 

The  plasma  frequency  cop  for  glass  is  small  compared  to  the  photon  frequency,  so  that  n  is  slightly  less 
than  unity.  Therefore,  x  rays  traveling  in  vacuum  or  air  can  be  totally  externally  reflected  from  smooth 
glass  surfaces.  The  maximum  angle  between  the  x-ray  beam  and  glass  surface  for  total  reflection,  0C,  is 
approximately 

6  c  =  ^.mrad,  (2) 

where  E  is  the  photon  energy  in  keV.  The  higher  the  energy  of  photon,  the  smaller  is  the  critical  angle. 
For  this  reason,  profile  errors  that  change  the  angle  of  incidence  are  more  detrimental  at  high  energies. 


1.D  Numerical  simulations 

To  evaluate  the  experimental  performance  of  capillaries  and  optics,  and  to  design  capillary  optics,  it  is 
necessary  to  be  able  to  predict  theoretical  behavior  for  complex  geometries.  Extensive  modeling  pro¬ 
grams  describing  the  propagation  of  x-rays  along  capillaries  and  optics  are  being  developed.  These 
modeling  programs  provide  essential  information  on  the  transmission  efficiency  and  divergence  of  cap¬ 
illary  optics,  and  are  being  used  to  simulate  the  performance  of  the  optics  by  obtaining  the  following 
outputs:  transmission  efficiency,  exit  divergence,  output  uniformity,  and  transmission  as  a  function  of 
entrance  angle. 


The  models  for  these  simulation  programs  are  based  on  Monte  Carlo  simulations  of  simple  geometrical 
optics.  The  simulations  randomly  generate  millions  of  x-ray  photons  from  a  particular  location  and 
solid  angle,  and  then  trace  them  through  the  optics  using  reflection  theory  and  classic  mechanics.  The 
necessary  confidence  in  the  simulations  is  being  built  by  comparing  the  results  of  the  modeling  program 
with  experimental  data. 


For  single  fibers,  the  computational  speed  is  greatly  enhanced  by  a  reduction  to  two  dimensions  by  pro¬ 
jecting  the  trajectory  onto  the  local  fiber  cross-section.8  Reflections  are  computed  from  standard  tables.9 
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1.E  Purpose 

The  purpose  of  this  project  is  to  contribute  to  an  ongoing  effort  to 
use  the  new  innovative  technologies  in  X-ray  optics  to  develop  a 
digital  mammographic  imaging  system  with  high  resolution  and 
contrast  by  extending  the  capabilities  of  the  computer  simulation. 

Enhanced  computer  simulations  are  critical  because  manufactur¬ 
ing  large  polycapillary  optics  is  difficult,  expensive  and  time- 
consuming,  and  simulation  is  a  good  design  tool  to  predict  the 
performance  of  the  optics  and  make  it  clear  how  all  the  different 
parameters  work  together  and  impact  the  resultant  transmission, 
resolution  and  contrast. 

1. F  Technical  objectives 

The  specific  aims  of  this  traineeship  are: 

•  to  simulate  all  trial  optic  setups,  to  compare  the  results  with  figure  3.  The  generation  of  a  ran- 
experimental  data  and  to  determine  the  best  design  for  build-  dom  point  (of  a  photon)  in  a  circle. 
ing  the  best  direct  digital  mammographic  system; 

•  to  further  the  author’s  knowledge  about  the  mammographic  systems  and  to  promote  a  future  career 
in  this  field. 

2.  Body:  Methods  and  Results 

2.A  Introduction  to  profile  and  surface  defect  corrections 

Real  glass  surfaces  are  not  perfect.  In  addition  to  absorption,  there  is  roughness,  which  has  short  wave¬ 
lengths,  of  the  order  of  microns  or  less,  waviness,  which  has  longer  wavelength,  and  bending,  which  has 
wavelength  longer  than  the  length  of  the  capillary.  Bending  is  described  by  the  radius  of  curvature  of 
the  fiber,  R  in  Table  1 .  Waviness  results  in  the  channel  wall  being  tilted  with  respect  to  the  axis  of  the 
fiber,  and  therefore  changes  the  angle  of  incidence.  It  is  described  in  terms  of  a  distribution  of  tilt  an¬ 
gles.  Surface  roughness  produces  both  diffuse  scattering,  for  which  the  reflected  angle  can  be  either 
larger  or  smaller  than  the  corresponding  incident  angle,  and  transmission  into  the  glass.  The  roughness 
correction10  used  in  previous  simulations11  is  described  in  terms  of  a  roughness  height,  z  and  a  correla¬ 
tion  length,  s. 

2.B  Current  progress 

To  provide  a  framework  for  assessing  the  results  of  the  first  year  of  the  project,  the  original  statement  of 
work  for  the  full  three-year  project  is  reproduced  below. 


Proposed  Statement  OfWork 


1)  Developing  and  testing  the  simulation  software,  (months  1-6) 

2)  Simulations  of  Post-patient  line  optics  matched  with  linear  array  detector,  (months  7-9) 

3)  Simulations  of  Post-Patient  Two  Dimensional  Optic  matched  with  Two  Dimensional  Detector, 
(month  10-12) 

4)  Simulations  of  Pre-patient  Line  Source  with  Post-Patient  Grid,  (months  13-14) 

5)  Simulations  of  Pre-patient  Line  Source  without  Post-Patient  Grid  (month  1 5). 
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6)  Analysis  of  the  results  from  the  tasks  2-5  and  selection  of  the  final  prototype  optics,  (months  16-18) 

7)  Design  of  final  prototype  optics,  (months  19-23) 

8)  Simulation  and  Comparison  with  the 
experimental  measurements  of  the  fi¬ 
nal  prototype  optics  and  analysis  of 
the  impact  of  non-ideal  performance 
on  mammographic  imaging,  (month 
24-35) 

9)  Final  report  of  the  project  (month  36). 


In  this,  the  first  year  of  this  project,  we 
developed  and  tested  a  more  physical 
model  and  compared  the  simulation  re¬ 
sults  for  a  variety  of  straight  capillaries  to 
measured  data  in  the  energy  range  from 
10  to  80  keV.  Using  this  analysis,  the  number  of  required  parameters  in  this  energy  range,  including  20 
keV  for  mammography,  is  reduced  from  four  to  two,  bending  curvature  and  waviness.  It  was  not  found 
to  be  necessary  to  include  roughness  effects. 


Figure  4.  Scheme  of  an  x-ray  beam  interactions  with  a 
randomly  tilted  surface. 


Although  It  took  almost  a  year  to  develop  this  unexpected  new  model  and  its  implementation,  it  was 
worthwhile  because  we  now  have  a  more  physical  simulation  model  and  it  will  increase  the  time  effi¬ 
ciency  to  do  tasks  2-8  in  the  future. 

The  following  sections  are  the  comparison  of  the  new  model  with  the  previous  model  and  the  compari¬ 
son  of  their  simulation  results  with  the  experimental  data. 

2.  C  Improvement  in  generating  x-ray  photons 

As  shown  in  Figure  3,  to  produce  a  random  point  P  in  a  circle,  we  can  randomly  generate  either  (r,  0)  or 
(x,  y).  Previous  simulations1 1  used  uniformly  generated  0  in  the  range  (0,  2n)  and  r  in  the  range  (0,  R), 
where  R  is  the  radius  of  the  x-ray  source,  capillary  channel  or  fiber  input  area.  This  work  uniformly 
generated  x  and  y  coordinates,  with  the  limitation  x2  +  y2  <  R2,  which  is  closer  to  an  isotropic  distribu¬ 
tion. 

2.D  Improvements  in  waviness  correction  model 

Waviness  of  a  capillary  surface  is  oscillations  with  wavelengths  shorter  than  the  capillary  length  and 
longer  than  the  wavelength  of  the  roughness.  The  waviness  effect  is  modeled  as  a  random  tilt  of  the 
glass  wall,  which  results  in  a  random  change  in  the  reflection  angle  on  each  bounce.  If  the  random  tilt 
angle  is  <|>  and  the  incident  angle  is  0,  then  the  reflection  angle  with  respect  to  the  nominal  surface  will 
be  0+2<|),  as  shown  in  Figure  4.  Previous  simulations  used  a  uniform  distribution  model;  this  work  uses  a 
tilt-corrected  normal  distribution  to  generate  <j). 

2.D.1  Uniform  distribution  model 

Previous  simulations,  referred  here  as  Model  Ml,  assumed  a  uniform  distribution  of  2<(>  from  -A0maxto 
A0max  if  0  ^  A0max,  or  from  -0  to  A0max  if  0  <  A0max.  The  maximum  random  tilt  angle  A0max  is  the  adjust¬ 
able  parameter  which  describes  the  amount  of  waviness  of  the  capillary. 
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2.D.2  Tilt-corrected  normal  distribution  model 


2.D.2.1  Normal  distribution 


For  this  work,  named  as  Model  M2,  it  is  assumed  that  these  tilt  angles,  <|>,  are  normally  distributed  in  the 
range  (-7t/2, rc/2).  with  the  mean  value  equal  to  zero.  For  high  quality  optics,  the  standard  deviation  of 
this  normal  distribution,  a,  is  much  smaller  than  the  critical  angle,  0C.  The  probability  distribution  of  tilt 
angles,  is 


G(4>)  = 


2<r2 


(3) 


2.D.2.2  X-ray  impact  tilt  correction 

In  this  work,  Model  M2,  consideration  was 
taken  of  the  fact  that  the  surface  tilt  angle  will 
affect  the  probability  of  x-ray  impact  on  that 
surface.  Taken  to  extremes,  a  surface  region 
perpendicular  to  the  beam  is  much  more  likely 
to  intercept  the  beam  than  a  surface  region 
parallel  to  the  beam.  Figure  5  displays  three 
surfaces,  OAi,  OA2  and  OA3,  with  different 
tilt  angles,  <J>i,  <j>2  and  fo,  respectively,  from  the 
nominal  surface  00’.  The  projections  onto  the 
nominal  surface  for  the  three  surfaces  are 
equal,  OAi*cos<|>i  =  OA2*cos<|>2  =  OA3*cos<|>3 

5 

The  probability  of  incidence  of  a  parallel  x- 
ray  beam  with  incident  angle  0  (with  respect 
to  the  nominal  surface  00’)  hitting  these 
tilted  surfaces  is  given  by  their  corresponding 
perpendicular  length  AjBj  (for  j  =  1, 2  or  3) 

P  oc  A.B.  =  OA,  ■  sin(i9 +<)>,)=  ^  ■  sin(0 + tf> ), 

111  1  COS0J  1 


Figure  5.  Three  surfaces,  OAj,  OA2  and  OA3,  with 
different  tilt  angles,  <f>i,  fa  and  fo,  respectively,  from 
the  nominal  surface  00  Here,  -6  <  <j>, 3  <0, 
OAj*cos  (/>  1  =  OA2*cos  (j> 2  —  OAfrcos  (f> 3  =  OO’. 


We  call  this  the  tilt-corrected  probability  distribution.  The  complete  description  of  this  distribution  is 


//(«>)  = 


F  sin(#  +  <t>) 
cos(  (j>) 

0, 


2 


(5) 


where  0  is  the  incident  angle,  <))  is  the  tilt  angle  and  F  is  a  normalization  constant. 


2.D.2.3  Tilt-corrected  normal  distribution 

Combining  the  normal  distribution  G(<|))  with  the  tilt  correction  H(<|))  gives  <|)  for  a  certain  incident  angle 
0  as  follows: 
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m= 


£  *+*> 


Ke^  • 


cos^ 


0, 


where  0  is  the  incident  angle,  §  is  the  tilt  angle  and  K  is  a  normalization  constant. 
Noting  that 

sin^+^)  =sin0+cos0-tan^ 
cos^ 


(6) 


(7) 


and  that  both  the  incident  angle,  0,  and  the  tilt  angle,  <|>,  are  very  much  less  than  40  mrad,  we  use  the  ap¬ 
proximation 


P(</>)  = 


It 


Ce  2<t  *(0  +  <j>),  -e<<j><4- 


0, 


Tt 


- <<t><-6 

2 


where  the  normalization  constant ,  C,  is 


C 


■Jin  9c t  +  2 <r 2 


(8) 


(9) 


to  speed  the  calculation.  P(<j>)  is  called  the  tilt-corrected  normal  distribution.  In  model  M2,  c  is  the  vari¬ 
able  that  describes  the  amount  of  waviness.  Two  examples  of  J(§)  and  P(<)>)  are  shown  in  Figure  6  and 
Figure  7. 


Figure  6.  The  probability  distribution,  J(<j>), 
and  its  approximation,  P (</>),  versus  tilt  angle 
<f>  when  incident  angle  9=  0. 0001  rad,  stan¬ 
dard  deviation  <7  =  0. 0004  rad.  <f>  is  in  rad. 


Figure  7.  J(</>)  and  P(<j>)  versus  tilt  angle  <f> 
when  incident  angle  6  =  0.009  rad,  standard 
deviation  <7  =  0. 0004  rad.  <f>  is  in  rad. 
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Fiber  Description 

Model  Ml 

Model  M2 

Fiber 

# 

Type 

Channel 

Size 

pm 

Open 

Area 

Length 

mm 

z 

nm 

s 

(im 

R 

m 

AOmax 

mrad 

R 

m 

a 

mrad 

A 

1 

0.5 

12 

65% 

105 

0.7 

6 

105 

0.4 

60 

0.225 

C 

3 

0.75 

22 

50% 

136 

0.5 

6 

125 

0.35 

225 

0.2 

D 

4 

4 

12 

55% 

130 

0.8 

6 

110 

0.285 

139 

0.125 

E 

5 

0.3 

55% 

105 

0.7 

6 

28 

0.2 

31 

0.09 

F 

4 

4 

12 

55% 

130 

0.8 

6 

90 

0.45 

90 

0.18 

Table  1  Parameters  for  best-fit  simulations.  R  is  the  bending  radius.  For  Ml,  A0max  is  the  amount 
of  waviness,  z  is  the  roughness  height  and  s  is  the  roughness  correlation  length.  For  M2,  <j  is  the 
standard  deviation  of  the  waviness. 


2.E  Results 

Transmission  of  a  number  of  different  fibers  measured  previously11  in  the  energy  region  10-80  keV  are 
shown  in  figures  8-18.  The  descriptions  of  these  polycapillary  fibers  and  their  simulation  parameters  are 
shown  in  Table  1. 
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Figure  8  Transmission  spectra 
for  fiber  C  simulated  with  dif¬ 
ferent  bending  curvature  alone 
using  models  Ml  and  M2, 
compared  with  experimental 
data.  This  figure  shows  that 
the  simulations  with  bending 
alone  can  not  fit  the  data,  and 
indicates  that  bending  is  not 
the  only  factor  to  cause  the 
transmission  drop. 


Figure  9  Simulations  of 
transmission  spectra  for  fiber 
C  with  only  waviness  com¬ 
pared  with  the  experimental 
data.  This  figure  shows  the 
effects  of  waviness.  It  also 
shows  simulations  using 
waviness  alone  do  not  fit  the 
data.  The  simulations  do  not 
include  the  roughness  or 
bending  correction. 


Figure  10  Simulations  of  trans¬ 
mission  spectra  with  waviness 
correction  using  different  numbers 
of  photons,  compared  with  the  ex¬ 
perimental  data  of  ’ fiber  C.  The 
figure  shows  that  the  simulation 
curves  using  10000  photons  are 
much  smoother  than  those  using 
1000  photons  because  of  the  re¬ 
duction  in  statistical  fluctuation. 
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The  experimental  data  and  simulation  results  from  the 
two  simulations  for  a  type  3  polycapillary,  Fiber  C,  are 
compared  first  to  demonstrate  the  fitting  processes,  to 
interpret  the  waviness  and  bending  effects,  and  to  un¬ 
derstand  the  final  results. 

2.E.1  Bending  effects 

A  slight  bending  can  dramatically  reduce  the  transmis¬ 
sion  of  high-energy  photons  because  of  the  small  criti¬ 
cal  angle  at  high  energies.  Figure  8  shows  a  compari¬ 
son  between  experimental  data  and  simulations  with 
different  bend  radii.  The  difference  in  the  bending  ef¬ 
fect  between  Ml  and  M2  is  due  to  the  different  algo¬ 
rithms  used  in  generating  x-ray  photons  as  discussed  in 
section  2.C.  Ml  assumes  more  photons  hit  the  center  of 
the  fiber  as  well  as  the  center  of  each  channel.  M2 
generates  photons  uniformly  across  the  channel,  fiber 
and  source. 

For  both  models,  a  bend  of  R  =  80  m  reduces  the 
transmission  too  much  at  high  energies.  Since  any 
waviness  or  roughness  correction,  like  the  bending  cor¬ 
rection,  will  reduce  the  transmission,  R  must  be  larger 
than  80  meters.  We  set  Rmjn  =  80  m,  which  is  useful  in 
determining  the  best-fit  parameters  for  the  model  M2. 


Figure  11.  Simulated  transmission  spectra 
using  model  M2  with  fixed  waviness  (<j=0.2 
mrad)  and  different  bending,  compared  with 
the  experimental  data.  The  simulation  curves 
shown  are  just  a  few  representative  cases  in 
the  bisection  process.  For  M2,  the  best-fit 
parameters  are  a  =0.2  mrad  and  R  =  225 
m;  for  Ml,  they  are  z  =  0.5  nm,  s  =  6  pm, 
AQmax  -  0.35  mrad  and  R  =  125  m. 


Figure  12.  Simulations  of  source  scan  curve  us¬ 
ing  model  M2  with  different  waviness  (a)  com¬ 
pared  with  the  experimental  data.  The  larger  the 
waviness,  the  narrower  the  simulated  scan  curve. 
The  simulation  curves  with  the  best  fitting  pa¬ 
rameters  of  Ml  and  M2  obtained  in  Figure  1 1 
are  also  shown  in  this  figure  for  reference. 


Source  positjon(mm) 

Figure  13  Simulations  of  source  scan  curves 
using  model  M2  with  no  waviness  and  differ¬ 
ent  bending  compared  with  the  experimental 
data  at  10  keV.  The  simulated  scan  curve 
does  not  become  narrower  even  when  the 
bending  radius  R  =  20  m.  The  word  “ideal” 
refers  to  a  simulation  with  R=ooand  no 
waviness. 
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Clearly,  bending  alone  cannot  explain  the  experimental 
data,  especially  at  low  energies. 

2.E.2  Waviness  effects 

To  explain  the  data  at  energies  from  20  to  60  keV,  an¬ 
other  correction  is  needed.  Roughness  does  not  have  a 
major  effect  in  this  energy  range11.  Waviness  is  ex¬ 
pected  to  be  the  major  factor  in  this  energy  range. 

Figure  9  shows  that  the  transmission  becomes  lower 
for  all  the  energies  and  the  high-energy  transmission 
drops  faster  with  increasing  a  gets  larger.  It  also  shows 
that  waviness  alone,  like  bending  alone,  can  not  ex¬ 
plain  the  experimental  data.  In  the  simulation,  a  =  0.3 
mrad  and  0.4  mrad  are  definitely  over-correcting  at  low 
energies,  so  we  can  set  0.3  mrad  as  the  maximum  of  c 
(cimax)  for  Fiber  C.  The  waviness  parameter  A0  for 
model  Ml  is  larger  than  the  parameter  a  for  model  2 
because  the  tilt  correction  in  model  M2  makes  it  more 
probable  that  the  x-ray  will  intersect  the  high  angle  sur¬ 
faces. 

The  fluctuation  of  the  simulation  curve  in  Figure  9 
when  a=0.25  mrad  is  mainly  due  to  the  relatively  small 
number  of  photons,  1000,  used  in  the  simulation.  Figure  10  shows  the  difference  when  the  number  of 
photons  changes  from  1000  to  10000. 

Using  a  bisection  algorithm  with  Rmjn  and  CTmax,  R  and  a  were  applied  together  and  compared  to  the  data 
to  get  the  best  fit  parameters,  as  shown  in  Figure  1 1 .  The  best  fit  parameters  for  fiber  C  using  model  M2 
are  R  =  225  m  and  cr  =  0.2  mrad. 

2.E.3  Verification  using  the  source  scan  curves 

Figure  12  shows  the  source  scan  curve,  that  is  transmission  as  a  function  of  lateral  source  displacement, 
and  simulation  results,  for  increasing  waviness  for  Fiber  C  at  10  keV.  The  simulated  source  scan  slowly 
becomes  narrower  when  a  gets  larger,  however  that  change  is  too  small  to  be  used  to  determine  a. 
Figure  13  displays  simulations  for  Fiber  C  at  10  keV  with  different  R  and  no  waviness.  Figure  14 
shows  the  effect  of  different  R  with  a  fixed  <j.  From  these  two  figures,  it  is  apparent  that  changing  R 
from  20  meters  to  infinity  does  not  have  any  visible  effect  on  the  fitting  curve.  Therefore,  we  can  not  use 
the  simulations  of  the  source  scan  curves  to  determine  either  a  or  R,  but  they  can  be  used  to  verify  a. 

For  Fiber  C,  as  shown  for  10  keV  in  Figure  12,  at  higher  energies  in  Figure  15,  and  as  a  function  of  en¬ 
ergy  in  Figure  1 1 ,  the  simulation  curve  of  a=0.2  mrad  fits  the  experiment  data  very  well  with  one  fewer 
fitting  parameter  than  the  best  fit  of  Ml . 


Source  position(mm) 

Figure  14  Simulations  of  source  scan  curves 
using  model  M2  with  fixed  waviness  (a  =  0.2 
mrad)  and  different  bending  compared  with 
the  experimental  data.  The  simulated  scan 
curve  does  not  apparently  change  even  when 
the  bending  radius  R  =  20  m.  With  and  without 
the  bending,  simulations  using  cr  =0.2  mrad fit 
the  experimental  data  quite  well. 
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20KeV  40KeV:  59KeV:  80KeV: 

*  data  *  data  *  data  *  data 

--■--Ml:  Best  fit  —"—Ml:  Best  fit  M 1 ;  Best  fit  Ml:  Best  fit 

•■■o---  M2:<y-0.2  “ -o -  M2:ct=0.2  --O-- M2:a=0.2  - -O- -  M2:ass0.2 

—A—  M2:a=0.2,R=225  •  A  •  •  M2:a=0.2,R=225  •  A  - M2:a=0.2,R=225  M2:a*0.2,R*225 


Source  position(mm)  Source  position(mm)  Source  position(mm)  Source  position(inm) 

Figure  15  Simulated  source  scan  curves  using  Ml  and  M2  compared  with  experimental  data  at  four 
different  photon  energies.  They  all  fit  the  data  quite  well.  For  M2,  this  figure  also  verifies  that  the 
bending  has  no  visible  effect  on  the  source  scan.  The  best-fit  parameters  for  Ml  are  R  =  125  m,  A0max 
=  0.35  mrad,  s  =  6  pm  and z  =  0.5  nm. 

2.E.4  Best-fits  of  all  the  fibers  and  analysis  of  the  results 


Using  the  algorithm  described  above,  the 
best-fit  parameters  of  M2  for  Fiber  A,  D 
and  E  were  obtained.  The  best-fit  parame¬ 
ters  for  both  Ml  and  M2  are  shown  in  Ta¬ 
ble  1.  The  best-fit  curves  of  Ml  and  M2 
along  with  the  experimental  data  are  shown 
in  Figures  16-18.  The  fibers  are  described 
in  Table  1. 

In  Figure  17,  the  transmission  for  fiber  A 
shows  a  rapid  drop  for  energies  above  30 
keV.  Although  fiber  D  has  lower  fractional 
open  area  than  fiber  A,  its  transmission  ex¬ 
ceeds  that  of  fiber  A  at  energies  above  30 
keV.  This  is  because  fiber  A  is  thin  (0.5 
mm  in  outer  diameter)  and  flexible,  there¬ 
fore  difficult  to  keep  straight  in  the  meas¬ 
urement  apparatus.  Model  M2,  which  can 
vary  only  bending  and  waviness,  requires  a 
much  sharper  bend  for  fiber  A  than  for  fi- 


Photon  energy  (KeV) 

Figure  16.  Simulations  of  transmission  spectra  of fiber 
C,  fiber  D  and fiber  E  with  their  best-fit  parameters 
compared  with  the  experimental  data.  All  parameters 
related  to  this  figure  are  listed  in  Table  1. 
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ber  D.  This  bending  is  more  significant  at  high  energies,  where  the  critical  angles  are  smaller. 

Transmission  curves  for  fiber  C  and  D  are  similar  in  shape  and  are  nearly  flat  up  to  60  keV,  as  shown  in 
figure  15.  However,  the  smaller  channel  size  for  fiber  D  results  in  the  simulation  being  less  sensitive  to 
bending.  Therefore  a  larger  bending  curvature  (smaller  radius)  is  required  for  fiber  D.  The  waviness 
correction  for  fiber  D  must  therefore  be  smaller  to  keep  the  transmission  about  the  same  as  for  fiber  C. 

Fiber  E  is  the  thinnest  fiber.  Table  1  shows  that  it  has  the  largest  bending  correction  because  of  its  flexi¬ 
bility.  Its  transmission  curve  is  flat  up  to  40  keV.  This  is  due  to  its  very  small  channel  size,  as  discussed 
in  the  comparison  of  fiber  C  and  fiber  D.  However,  if  the  channel  size  is  too  small,  as  in  fiber  E,  it  also 
results  in  more  reflections  being  needed  for  a  photon  to  traverse  the  fiber  and  may  have  introduced  other 
defects  such  as  blocked  channels.  This  is  why  the  transmission  is  only  40%  for  the  energies  below  40 
keV  although  the  open  area  is  around  55%. 


The  high  transmission  and  the  simulation  results 
good.  The  bending  radius  is  above  130  meters 
for  type  3  and  type  4  capillaries.  It  is  hopeful 
that  we  can  further  improve  the  high-energy 
transmission  performance  of  polycapillary  fi¬ 
bers  by  decreasing  to  an  optimum  channel  size, 
and  making  them  more  rigid. 

Because  Model  M2  is  physically  plausible, 
gives  reasonable  results,  and  matches  the  ex¬ 
perimental  data  very  well,  as  shown  in  figures 
16  and  17,  we  can  use  it  to  simulate  other  poly- 


Figure  18.  Simulations  of  transmission  spectra 
of  fiber  D  (unexposed)  and fiber  F  (exposed) 
with  their  best-fit  parameters  compared  with 
the  experimental  data.  Layer  =  35  means  that 
there  is  a  35  micron  glass  layer  added.  All  pa¬ 
rameters  related  to  this  figure  are  listed  in  Ta¬ 
ble  1  except  the  glass  layer. 


show  that  the  quality  of  the  capillary  fibers  is  quite 


Figure  17.  Simulations  of  transmission  spectra  of 
fiber  A  and fiber  D  with  their  best-fit  parameters 
compared  with  the  experimental  data.  All  parame¬ 
ters  related  to  this  figure  are  listed  in  Table  1. 

capillary  optics.  Compared  to  Ml,  which  uses  four 
parameters,  M2  uses  only  two  parameters,  which 
makes  simulating  the  experimental  results  more 
straightforward,  easier  and  timesaving.  In  addition, 
roughness  does  not  have  much  effect  in  the  energy 
range  10-80  keV,  so  M2  gives  us  more  confidence 
with  fewer  parameters. 

2.E.5  An  analysis  of  a  radiation-damaged  fiber 

The  measured  transmission  spectra  and  source  scans 
provide  sensitive  tools  for  analyzing  capillary  quality, 
This  method  was  used  to  analyze  the  effect  of  radia¬ 
tion  damage  on  fiber  F  that  has  the  same  type  and 
same  length  as  fiber  D.  Radiation  damage  has  been 
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observed  after  a  very  large  dose  of  x-ray  exposure.12  It  was  found  that  prolonged  exposure  to  synchro¬ 
tron  white  beam  radiation  could  cause  a  measurable  bend  for  a  thin  fiber.  This  bending  can  be  reduced 
by  annealing  and  by  holding  the  fiber  rigid  during  exposure.  Fiber  F  had  1.8  MJ/cm1 2 3 4 5 6 7 8 9 10 11  exposure  before  the 
current  measurements.  There  is  no  visually  observable  bending  or  other  visible  damage. 

The  transmission  spectrum  and  its  simulation  for  fiber  F  are  compared  with  an  unexposed  capillary  of 
the  same  type,  fiber  D,  in  Figure  16.  The  measured  data  shows  a  significant  transmission  drop  above  30 
keV  and  below  10  keV.  The  simulation  indicates  that  the  exposed  capillary  has  more  bending  and 
waviness.  The  bending  radius  was  changed  from  139  m  to  90  m  and  the  waviness  a  changed  from  0.125 
mrad  to  0.18  mrad.  Although  these  profile  changes  have  a  significant  effect  the  transmission,  they  are 
still  very  small  and  not  visible.  The  transmission  drop  below  20  keV  can  be  explained  by  adding  35  // m 
x-ray-absorbing  glass  layer  as  shown  in  Figure  18.13 

3.  Conclusions 

A  new,  more  physical  model  for  describing  profile  defects  in  polycapillary  x-ray  optics  at  high  energies 
for  mammography  has  been  developed.  The  simulation  model  M2,  a  geometric  simulation  with  wavi¬ 
ness  and  bending  corrections,  is  physically  plausible  and  works  quite  well  in  explaining  the  transmission 
spectra.  It  uses  only  two  parameters,  which  makes  simulating  the  experimental  results  and  designing 
lenses  straightforward,  easy  and  timesaving.  According  to  the  simulations,  waviness  and  bending  in  the 
capillary  channel  profiles  can  be  particularly  harmful  for  high-energy  photons  because  of  the  smaller 
critical  angle  at  high  energies.  Roughness  is  not  important  in  the  simulation.  It  is  expected  that  increas¬ 
ing  the  rigidity  and  straightness  of  the  fiber  and  reducing  the  channel  size  to  an  optimum  value  can  fur¬ 
ther  increase  the  transmission.  The  spectra  of  transmission  versus  energy  have  been  found  to  be  a  very 
sensitive  tool  in  capillary  quality  analysis,  and  transmission  versus  source  scan  curve  measurements  can 
also  provide  significant  information. 
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